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O ! Abstract 

o 

In the present letter, a particular form of Slavnov- Taylor identities for the Curci-Ferrari model 
. is deduced. This model consist of Yang-Mills theory in a particular non-linear covariant gauge, 

supplemented with mass terms for gluons and ghosts. It can be used as a regularization for the 
Yang-Mills theory preserving simple Slavnov- Taylor identities. Employing these identities two non- 

■ renormalization theorems are proved that reduce the number of independent renormalization factors 
from five to three. These new relations are verified by comparing to the already known three-loops 
renormalization factors. These relations include, as a particular case, the corresponding known 

■ identities in Yang-Mills theory in Landau gauge. 
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1 Introduction 



In last years important progress has been made in the understanding of the non- 
perturbative sector of Yang-Mills theory and QCD. The main reason for that comes 
from extensive Monte-Carlo simulations that have made possible the calculation of dif- 
ferent non-perturbative quantities. Lattice simulations have so opened the possibility to 
test non-perturbative analytical approximations. Probably the most important results of 
this kind are the Yang-Mills correlators, mainly using Schwinger- Dyson [1-14] and Non- 
Perturbative Renormalization Group [15, 16] equations and which have been compared 
quite successfully with lattice results [17-23]. These gauge-dependent objects have been 
calculated in a variety of gauges (see, for example, [24-27]) but most results have been 
obtained in Landau gauge. 

These results indicate that analytical calculations with relatively simple approxima- 
tions have a chance to be controlled not only in the ultraviolet sector of the theory but 
also in the infrared, at least in Landau gauge. Contrarily to what happens in perturbation 
theory, the strong coupling does not seem to diverge in the infrared at a finite momentum 
scale. Moreover, even if results depend on the scheme and on the details of the approx- 
imations used, different analyses seem to agree on the fact that many naturally defined 
Yang-Mills couplings remain bounded in the infrared. In the context of Landau gauge, 
the ghost-gluon coupling becomes particularly simple to analyze due to the well-known 
non-renormalization theorem for this specific coupling [28] . Many lattice results even seem 
to indicate that this coupling goes to zero in the infrared [22, 23]. On the same lines, gluon 
propagator do not have a massless behaviour as obtained in perturbation theory but a 
massive one (Schwinger-Dyson indicates generally an infinite mass [1,2,4,7,9,10] and 
lattices a finite one [17, 18, 22, 23]). 

This massive behaviour in the gluonic sector, with a frozen coupling in the infrared, 
reduces considerably infrared instabilities of Yang-Mills theory. However, the ghost sector 
remains massless in Landau gauge at a non-perturbative level and analytical or semi- 
analytical treatments must face typical difficulties of critical phenomena. In most cases 
where quantitative results for critical systems have been achieved some sort of renormal- 
ization group procedure has been necessary. The reason is that renormalization group 
equations only involve modes in a relatively small shell of momenta. This is particularly 
true in the Wilsonian version of the renormalization group [29-36] but in some sense is 
the characteristic property of different kinds of renormalization group equations. This 
must be contrasted with Schwinger-Dyson equations where all momenta are coupled in 
a significant way. One manifestation of this is that these equations require an explicit 
renormalization. To implement this renormalization process in a non-perturbative way is 
far from trivial. This is a serious motivation in order to look at the construction of renor- 
malization group equations in Yang-Mills theory describing in a complete way correlation 
functions. 

The ideal tool in this sense would be to have a Wilsonian flow equation for Yang- 
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Mills theory. However, the standard procedure [37-44] does not preserve simple Slavnov- 
Taylor identities [28, 45] along the flow and an alternative formulation that preserve gauge 
symmetry has turned to be extremely involved [46-52]. Another renormalization group 
given by Callan-Symanzik equation [53, 54] describes the flow of a field theory in term of 
a renormalized mass parameter that regulates completelly the theory in the infrared. It 
has poorer non-perturbative properties than Wilsonian flows, but is much simpler. In the 
Yang-Mills case, this requires a deformation of the theory in such a way that gluons and 
ghost acquire a mass, preserving renormalizability and reasonably simple Slavnov- Taylor 
identities. The simplest way to do that is offered by the Curci- Ferrari model [55,56]. It 
corresponds to the Yang-Mills theory in a certain family of gauge fixings supplemented 
with a particular mass terms for gluons and ghosts. If the gauge and mass parameters 
go to zero, Curci-Ferrari model reduces to the important case of Yang-Mills theory in 
Landau gauge. It may be useful to remember here that only the theory with zero mass 
parameter is physically acceptable. In fact, it has been proved that only in this case the 
model verifies perturbative unitarity and gauge parameter independence (see for example 
[56-58]). This, however, does not excludes the use of this model as an intermediate step 
regularization if at the end of the calculation the mass parameter is taken to be zero. 

There is another motivation to consider the Curci-Ferrari model. The non-perturbative 
definition of a BRST symmetry [59-61] has been elusive. In particular, Neuberger has 
constructed a lattice version of the BRST symmetry [62] but, he also showed [63] that 
averages of gauge invariant operators take the form in any finite lattice of a zero over zero 
expression. This difficulty is closely related to the Gribov problem [64] of non-perturbative 
gauge fixing in non abelian gauge theories. However, recently it has been proved [65] 
that the Curci-Ferrari model regularize this zero, gauge invariant correlators acquiring a 
definite expression for any non- zero mass. This allows for a non-perturbative construction 
of the BRST-invariant model on a lattice [66]. This is radically different to usual gauge 
theories where Neuberger's zero forbids such a construction [63]. In this way, Yang-Mills 
theory can be seen as the zero limit of the Curci-Ferrari model whose BRST symmetry 
has a non-perturbative meaning. Moreover, the corresponding limit can be seen as the 
result of the renormalization group flow. 

Finally, the Landau-gauge has more symmetries than Curci-Ferrari gauge and so is 
fixed under renormalization. However, if one takes a zero gauge parameter along the flow, 
ghost remain massless and Neuberger's zero is not regularized. Fortunately, the Landau 
gauge (where most of the gauge-fixed Monte-Carlo simulations have been performed) 
seems to be an attractive renormalization group fixed point (see, for example, [43]). If so, 
the zero mass limit automatically pushes the theory to the Landau gauge in the infrared. 

For all these reasons, the analysis of Curci-Ferrari model is interesting. However, even 
if multiplicatively renormalizable, Curci-Ferrari model seemed to require five renormal- 
ization factors [57]. In fact, a claim of a reduced number of renormalization factors [67], 
turned out to be incorrect [57]. These renormalization factors can be taken to be those 
corresponding to the renormalization of gluon and ghost fields, gauge parameter, coupling 
constant and mass and have been calculated up to three loops. This is to be compared 
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to usual linear covariant gauges where (without masses) three renormalization factors 
are necessary or to the Landau gauge where only two are required. This complicates 
considerably any analytical or semi-analytical analysis. In this letter, a perturbative anal- 
ysis of the renormalization procedure of the Curci-Ferrari model is performed showing 
that, in fact, only three renormalization factors are required. Unexpectedly, some sym- 
metries present in the model had not been fully taken into account up to now and, in 
fact, the mass and coupling renormalization can be expressed in terms of the other three 
renormalization factors. These relations are generalizations of the already known Lan- 
dau gauge non-renormalizations theorems for these quantities [28,68]. This result is only 
a consequence of the high level of symmetry of the model that may be very useful in 
non-perturbative approximations. 

The outline of the letter is the following: First, the Curci-Ferrari model and its symme- 
tries are reviewed. Second, the corresponding Slavnov- Taylor identities and the associated 
consequences for the renormalization factors are deduced. Finally some conclusions and 
perspectives are presented. 



2 The Curci-Ferrari model 



The Curci Ferrari model is a modification of Yang-Mills theory that includes masses for 
gluons and ghost preserving a sort of BRST symmetry. It has also many other symmetries, 
one of them being an anti-BRST one. Here the model without matter fields in a 4- 
dimensional euclidean space-time is considered, but it is straightforward to see that the 
results to be proved are true also including matter fields or considering a Minkowski 
space-time. 

The variations under BRST and anti-BRST of different fields are: 

sc a = -ff abc c b c c , sc a = -^dpA" - ff abc c b c c , (1) 

sc a = ^Al - ff abc c b c c , sc a = -ff abc c b c c . 

Here, g is the bare gauge coupling and £ is the bare gauge parameter of the model. A^ 
is the gauge field, c and c are ghost and antighosts respectively. D^c a denotes the ghost 
covariant derivative d^c a + g f abc A b c c and similarly for antighost. 

The most general renormalizable Lagrangian with these symmetries and the space 
time euclidean symmetry, global color invariance and the conservation of the number of 
ghost take the form: 

£ = £yM + £-GF + C>m- (2) 
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Here C Y m is the Yang-Mills lagrangian: 



£ym - -^Fj^F^, (3) 

where F® v = d^A^ — d u A a ^ + g f abc A b A^. Cgf is the gauge fixing term, that includes a 
ghost sector. It takes the form: 

Cgf = \d^D^c a + l -D^c a + ^(^) 2 - £ f f abc f ade c b c c c d c\ (4) 
Finally the mass term C m takes the form: 

c m = ml (±a;a; + £ ev) , (5) 

where m is the bare mass. 

It is important to observe that the action has also other linearly realized symmetries 
[69]. First, there is the discrete symmetry where c a — > c a and c a — > — c a . There are also 
two continuous symmetries. The first one has the infinitesimal form c a — > c a + ec a . The 
second, independent, symmetry correspond to c a — > c a + ec a . These symmetries, as BRST 
and anti-BRST do not present anomalies if the corresponding massless model do not. This 
may be seen from the fact that there is an explicit lattice regularization which respect all 
these symmetries [66]. 

BRST and anti-BRST transformations are not nilpotent in the Curci-Ferrari model 
even after imposing equations of motion. However, one can deduce closed Slavnov- Taylor 
identities by exploiting equations of motion for ghost and anti-ghost fields. Consider the 
s 2 and s 2 operations: 

s 2 Al = s 2 A* = s 2 c a = s 2 c a = 0, 

s 2 c a = -±-d,(D,c a ) + |/^c 6 ( - - ^f cde c d c e ^J + g Af^f^c d c e c c , 

Before deducing Slavnov- Taylor identities for these symmetries it is useful to consider also 
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Field/Source 


A 


c 


c 


K 


K 


L 


L 


M 


Dimension 


1 


1 


1 


2 


2 


2 


2 


2 


Ghost number 





1 


-1 


-1 


1 


-2 


2 






Table 1 

Canonical dimensionalities and ghost number of different fields and sources. 
'mixed' variations for different fields: 



ssA a a 



ssc 



ssc a 



ssc 



ssc a 



£>„( - jd v A a v - |/ abc c b c c ) + g f abc (D^)c c , 



ssA a u , 



90 f cde-d c e 



= 9of ahc c b (-U l/ Al-^f 

= -Ld^D^C a - OOf-abcf-cde-b^e 

£o 4 

,2 



SO , a 6e 

2 J 



6 2 J 



= ~d,D.C a + VAfabcfbde-d-e c + 90 f abc-bf _ }_Q _ SOyafegd e\ 



6 

So/ 



SO rcde-d^e 

2 ; 



(7) 



In order to deduce Slavnov- Taylor identities for these symmetries, it is necessary to intro- 
duce sources for the variations of different fields under BRST and anti-BRST symmetries. 
Because the operators are not nilpotent, one should include also sources for the s 2 and s 2 
variations. Moreover, the introduction of sources for the variations of fundamental fields 
under ss, should be necessary a priori (see, for exemple, [70]) because variations under 
anti-BRST transformations are not BRST invariant as just seen. However, as will be seen 
below, almost all these double variations can be expressed in terms of fundamental fields 
and simple variations by exploiting ghost and antighost equations of motion. The only 
variation that can not be obtained in this way is ssA®. Even this particular variation will 
not be needed in the present letter where Slavnov- Taylor identities are only considered at 
K = 0. Let's consider the generating functional: 

exp(H/[J, X , X , K, L,K,L]) = J D(A, c, c) exp j J d 4 x( - £ + Jffi + X a c a + X a c a 



+ K*sA* + L a sc a + K^sAl + L a sc a + M a (-g /2f abc c b c c ) 



- —M a M a + —L a L a + KX 
2£o Co 



mm 



(8) 



Notice that sources for composite operators are coupled non-linearly. As will be seen bellow 
this is necessary in order to preserve multiplicative renormalizability. Table 1 contains 
dimensionalities and ghost numbers for different fields and sources. As usual, the Slavnov- 
Taylor identity associated to BRST symmetry is obtained by performing the change of 
variables in the functional integral tp — > tp + Sstp [tp being all the fundamental fields and 
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5 a constant grassmanian parameter). One obtains for K = 0: 

J d A x(j; S Al - x a sc a - X a sc a + L a ssc a - |m" f abc s(c b c c ) 



K=0 



0. (9) 



In order to close the equation we need expressions for (sc a ), (ssc a ), and f abc (s(c b c c )) . To 
obtain them one can use the equations of motion for ghost and anti-ghost fields. The ghost 
equation of motion is obtained by performing in the generating functional the infinitesimal 
change of variables c a (x) — > c a (x) + e a (x) giving: 

-a i 90 rabc-b]\/rc , „ rabcfyrc r-> jsa I ^nA^« ^ A n tP- i ^ „2 rabc rcde-b-d „e\ < ^,2, 



X a +^r c c b M c +g f abc c b L c -D„KZ = {--D^-^ 

(10) 

Here the same notation has been used for the average of fundamental fields and for the 
fields themselves. Along the same lines one obtains: 

^ ' (11) 

This implies the following relations: 

U-s 2 c a ) = X a + jf abc c b M c + g f abc c b L c - D^K; - m 2 ^c a , 

Us 2 c a ) = X a + jf abc c b M c + g f abc c b L c - D,R; + m 2 ^c a , 
-U(srs + r ss )c a ) = x a + ^f abc c b M c + g f abc c b L c - D^K; - m^ c a , 

Ussc a ) = 2( - d^sAD -x a - jf abc c b M c - g f abc c b L c + D,K; + m^ c a ) , 
-U(ss + ss)c a ) = X a + ^f abc c b M c + g f abc c b L c - D,R; + m 2 ^ c a , 

Ussc a ) = 2(d»{sAD - X a - 9 jf abc c b M c - g f abc c b L c + D,R; - m^ c a ). 

(12) 

Finally, the last remaining variation can be obtained noting that: 

s( - f = (s(^d,Al - f /<^cV - 1^)) = (*V - ^sAfi). (13) 

Please observe an important non-trivial fact, all double variations (aside ssA 1 ^) have been 
expressed in terms of operators whose sources are included in eq. (8). 

Before considering renormalization properties of the model, it is important to realize 
that under the infinitesimal linear symmetry 

c a ^c a + ec\ (14) 
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composite operators mix in the following way: 
sAl -> sAl + esA" 

fl 

sc a -> sc a - eg f abc c b c c 

S (-f / abc c 6 c c ) - s ( - f / abc c 6 c c ) - e| / abc c b c c = 6sc a (15) 

all the other being invariant. This is equivalent to saying that sources for composite 
operators vary linearly in the following way: 

k; -> k; + eR; 

L a -> L a + eM a 

M a -*M a + 2eL a (16) 
On the same lines, one deduce that under the infinitesimal symmetry 

c a -> c a + ec a , (17) 
sources for composite operators vary linearly also: 

- + ex; 

L a -> L a + eM a 

M a -> M a + 2el a (18) 
Please observe that terms non-linear in sources are invariant under these transformations. 



3 Non-renormalization theorems 



In this section, the renormalization properties of the CF model implied by ST identities 
previously presented is considered. By substituting expressions (12) and (13) in (9) and 
expressing the results in terms of the effective action one arrives to: 



d 4 x< 



5T 5T 

6A*6K* 



5T 5T 



5T 



5T 



i ~~ ~~ _| f) A a t 

5c a 5L a 5c a \ £o 11 » 5M a 




-g f abc c b L c 



D K a 



K=0 



0. 



(19) 



The perturbative renormalizatibilty of CF model has been proven by considering 
five renormalization factors [57]. Here this proof is revisited by exploiting this particular 
form of the ST identity. The non-trivial use of equations of motion presented before will 
imply that only three renormalization factors are required. To prove that, I follow the 
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standard procedure [71,72] of considering terms that can diverge by power counting and 
constraining them iteratively. In a loop expansion, suppose that all divergences have been 
renormalized at order n. Diverging terms that appear in the order n + 1 in the effective 
action, have couplings with positive or zero dimension. Let's call them Ar^, and take 
an infinitesimal constant, e. If one calls r<n v = S + eATdiv, then, in dimensions lower than 
four, the most general form for this functional at order n + 1, takes the form: 

T dtv [A, c, c, K, K, L, L, M] = T[A, c, c] - J d d xi^l a L a - ^M a M a + 1 K;K« 

+ K^sAl + K*sA° + L a ~sc a + L a lc a + M a ^ - + ~sc a 

(20) 

where, for the moment, sA®, sA®, sc a , sc a and sc a denote arbitrary operators of dimension 
two with the same transformations under linear symmetries as the corresponding bare 
expressions. Here the coefficients of L a L a and M a M a are related by the linear symmetry 
(14). 

Substituting these expression into the Slavnov- Taylor identity (19), one can analyse 
first the particular case corresponding to L = M = 0. This gives the three relations: 

?Al = 0, (21) 
§V = 0, (22) 
sf = 0. (23) 

(see, for exemple, [71] for details of an analogous deduction in a similar situation). Now, 
the most general operators of dimension two, respecting Lorentz invariance, global color 
invariance, symmetries (14) and (17) and nilpotecy (21,22) are 

~sAl = Zd^ + gr hc A\c c , 
U% = Zd,c a + gf abc Alc c , 

~sc a = - 9 -j ahc c h c c , 
2 

lc a = -^f abc c b c c , 
2 

sc a = ]rd,Al - §r bc cV. (24) 
The equation linear inMatL = Z = fT = becomes then, 

^ + (25) 
Then, imposing the two identities (25,23) one deduces that the most general expression 
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for f takes the form f = T m + T GF + T Y m, where: 



f j4 m / 1 

= l dx z^Tz 



T GF = / d 4 J —d u ^D u c a + —D^d^c" + — *— (d„A") 2 - ^ °^ - f abc f ade c b c c c d c e 1, 



with 



1 



dM— (26) 





= d,c a + 


90 fabc \b „c 






= d,c a + 


90 rabc /ib -c 

Z 9 VZa~ J " ' 




pa 
[ii/ 


= d,A a v - 


£gV ^ A 


(27) 



These factors can be absorbed in the following renormalizations: 



Aa _ 


\f^A A R,^ 


c a = 


\fZ~ c C R 




c a = ■ 




6> = 




(28) 


m\ = 


Z m TP?fti 


9o = 


Zg9R. 





Previous parameters (20,24) can be expressed in terms of bare ones and renormalization 
factors (28) by 

Zm = f = zj^do = z^kl9R, 

where the relation to renormalized parameters has also been given. One concludes that 
renormalization factors are not independent. They must satisfy the non-renormalization 
theorem 

Z m = y±~. (30) 

Renormalization factors for the Curci- Ferrari model have been calculated up to three loops 
in the MS scheme by Gracey [73] and one can verify explicitly that these renormalization 
factors verify (30). In fact, in the recent article [74] this relation is observed in the three 
loops perturbative series but an all order argument and the associated symmetry is absent. 

Now, one can consider linear terms in L in the identity (19). Terms proportional to 
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L a L b c c give another non-renormalization theorem: 




(31) 



Notice that the Landau gauge is the intersection of usual covariant linear gauges and the 
Curci- Ferrari gauge. Then, in the Landau gauge limit £—►(), the covariant linear-gauge 
non-renormalization theorem of the gauge parameter is valid Zg ~ 1. Relations (30) and 
(31) then become usual Landau gauge's non-renormalization theorems [28,68]. 



4 Conclusions and perspectives 

In this paper some particular form of Slavnov- Taylor identities valid for the Curci- 
Ferrari model have been deduced. They allow to deduce two non renormalization theorems 
that reduce the number of independent renormalization factors in this model from five to 
three. These relations have been tested in the already known MS three loops calculation 
of these renormalization factors. One of the two identities implies that the mass renormal- 
ization is fixed by the renormalization of dimensionless coupling in d — 4. This simplifies 
considerably the corresponding Callan-Symanzik equation, whose analysis is actually in 
progress [75]. 

Even if in this paper this particular Slavnov- Taylor identity has only been exploited in 
order to study renormalization properties of the Curci-Ferrari model, it can also be used 
to simplify the analysis of finite parts of vertices of the model, mainly in order to develop 
non-perturbative renormalizations schemes in Schwinger-Dyson or similar contexts. 
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